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ABSTRACT 

We present results of high-cadence monitoring of the optical light curve of the nearby. Type la SN 2014J in M82 using the 2.3m 
Aristarchus telescope. B and V-band photometry on days 15-18 after t,„ndS)> obtained with a cadence of 2 min per band, reveals 
evidence for rapid variability at the 0.02-0.05 mag level on timescales of 15-60 min on all four nights, taking the red noise estimation 
at face value. The decline slope was measured to be steeper in the S-band than in R-band, and to steadily decrease in both bands 
from 0.15 mag day^' (night 1) to 0.04 mag day^' (night 4) in V and from 0.19 mag day^' (night 1) to 0.06 mag day^' (night 4) in B, 
corresponding to the onset of the secondary maximum. We propose that rapid variability could be due to one or a combination of the 
following scenarios: the dumpiness of the ejecta, their interaction with circumstellar material, the asymmetry of the explosion, or the 
mechanism causing the secondary maximum in the near-infrared light curve. We encourage the community to undertake high-cadence 
monitoring of future, nearby and bright supemovae to investigate the intraday behavior of their light curves. 
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1. Introduction 


Nearby supernovae offer the opportunity to explore the short- 
timescale and low-amplitude variability properties of their light 
curves. Even though several nearby, bright supernova^ (SNe) 
have been discovered recently (e.g. the Type la SN 201 Ife and 
SN 2013aa, the Type IIP SN 2013ej; |Nugent et al.|201 l||Waagen' 
2013[ |Richmond|2014| respectively), no high-cadence variabil¬ 
ity search has been undertaken so far to explore variability on 
timescales of minutes or hours. Typical photometric monitoring 
of SNe consists of a single observation per night or every few 
nights in several hlters, therefore, the intraday behavior of the 
light curves of SNe remains uncharted territory]^ 

To rectify the situation, we performed high-cadence photom¬ 
etry of the nearby Type la supernova (SN la) SN 2014J (see e.g. 
Eoley et al.|2014 and references therein), which was discovered 
in the galaxy M82 in January 2014 (|Eossey et al.|2014| l. Reach¬ 
ing Vmax - 10.61 + 0.05 mag ( |Eoley et al.|2014 l, it was ideally 
suited for a variability search, which we conducted in Eebm- 
ary 2014 with the 2.3m Aristarchos telescope. Just before the 
submission of this work, [Srverd et al.] ( |2015) l reported relatively 
high-cadence photometry of SN 2014J with the 4.2cm Kilode- 
gree Extremely Little Telescope North (KELT-N), thereby plac¬ 
ing a 4.5% (3cr) upper limit on short timescale variations. This 


* Based on observations made with the 2.3m Aristarchos telescope, 
Helmos Observatory, Greece, which is operated by the Institute for As¬ 
tronomy, Astrophysics, Space Applications and Remote Sensing, Na¬ 
tional Observatory of Athens, Greece. 

** Table 2 is only available in electronic form via 
http://www.edpsciences.org 

' http://www.rochesterastronomy.org/s nimages/snmag.html 
^ While this paper was under review. Oiling et al. (2015 i reported 
serendipitous, high-cadence monitoring (every 30 min) of three type 
la supemovae by the Kepler Mission. 


paper reports the results of our monitoring survey: the observa¬ 
tions and data reduction are presented in Secion 2, the analysis 
in Section 3, the results in Section 4, and the discussion and con¬ 
clusions in Section 5. 
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Eig. 1; V-band image of SN 2014J in M82, obtained on 2014 
Eebruary 16 with the Aristarchos telescope. The positions of the 
5 comparison stars (S1-S5) used in the analysis are labelled. 


Article number, page 1 of|^ 





















A&A proofs: manuscript no. ms 



CO 



Fig. 2; Calibrated F-band (circles) and B-band (triangles) differential light curves of SN 2014J, compared with Star 1, as a function 
of the time in days measured since lmax{B) - 2456690.5 (JD; Foley et al. 2014 1 . Each panel corresponds to one night of observations. 
The B-band curves are offset by 0.95 mag in the upper panels and 1.05 mag in the lower panels for clarity. The black lines represent 
least-square fits to the photometry. 


2. Observations and Data Reduction 

SN 2014J was monitored for variability with the LN CCD cam¬ 
era on the 2.3m Aristarchos telescope for 4 consecutive nights: 
2014 February 16-19 (nights 1^ or Nl^), corresponding to 
days 14.8-18.2 after W,,(B) = 2456690.5 + 0.2 (JD;|Foley et al. 
2^ . The 1024x1024 SITeAB back-illuminated CCD has pix¬ 
els that are 24 pm in size, which map to 0.28" pixel ' on the fo¬ 
cal plane, making each image 4.8' on the side. Figure[2provides 
a finder chart, labelling the position of the SN and the compari¬ 
son stars (S1-S5) used in the analysis below. 

SN 2014J was observed for about 8 hrs per night to ob¬ 
tain differential photometry. Sequences of alternating V-band 
and B-band images were obtained with 5 sec and 20 sec expo¬ 
sure times, respectively, yielding 400-500 images per band each 
night and a cadence of 2 min per band. Small gaps in the data 
are mainly due to recentering performed in between typical se¬ 
quences of 100 images, as the absence of guiding caused the 
target to drift during the sequences. 

Initial reduction of the images was performed using standard 
routines in the IRAlj^ ccBproc package (i.e. bias level subtrac¬ 
tion, flat field division). The images were aligned with the IRAF 


^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As¬ 
tronomy (AURA) under cooperative agreement with the National Sci¬ 
ence Foundation. 
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imalign task, while the seeing was calculated using the IRAF 
task psfrneasure. The statistical mode values for the seeing in 
the V-band images each night were 1.2", 1.4", 1.2", 1.2", re¬ 
spectively, and 1.2", 1.6", 1.05", 1.25" in the B-band. Typical cr 
values were 0.55" for N2 and 0.2 - 0.3" for the other nights. 
Note, on N1/N3, the S/N of the maximum pixel value of SI in B 
(190/185) exceeded the mean S/N of the B-flat (176/172). 


3. Analysis 


Aperture photometry was extracted from the supernova and the 
5 brightest comparison stars available in the field, using a 5 pixel 
or 1.4" radius with the IRAF apphot package and a 5 pixel wide 
sky annulus 15 pixels from the source. Typical photometric er¬ 
rors are 1 mmag for the SN and SI, 0.01 mag for S2, and 0.01- 
0.02 mag for S3-S5. 

The instrumental magnitudes for the SN and S1 were trans¬ 
formed to the standard system by computing zeropoint offsets 
and color terms, using B and V magnitudes for SN 2014J on N1 
from Foley et al. |2014 1 and for S1 from the APASS Catalogu^ 


The resulting transformation equations are: 


y = v- 1.320-H 0.056-(B-y) (1) 

* http://www.aavso.org/apass/ 
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B^b- 2.609 + 0.005 -{B-V) 


( 2 ) 


where b and v are the instrumental magnitudes. 

Differential light curves were created by taking the difference 
between the instrumental magnitude of the supernova and each 
of the 5 comparison stars. SI yields the smallest scatter, with an 
error in each point of 1.4 mmag. However, due to its brightness, 
it was saturated or exceeded the non-linearity limit of the CCD 
(~ 45000 ADU) in the best seeing Z?-band images, creating gaps 
in the differential curves with respect to S1. Furthermore, the ab¬ 
sence of guiding caused SI to intercept a bad column, creating 
spikes in the light curves in both bands, which were selected vi¬ 
sually and removed. Figurej^presents the differential, calibrated, 
light curves with respect to SI, which show the expected decline 
of the SN and onset of the secondary maximum discussed below. 

Figures [3]-^ show the differential light curves with respect 
to all comparison stars available each night. We used a detrend¬ 
ing algorithm to reconstruct the light curves, rule out intrinsic 
variability in S1 (as no variability information for S1 is available 
in the literature), and correct for systematic effects. We applied 
the Trend Filtering Algorithm (TFA; Kovacs et al.|2005), as im¬ 


plemented in the VARTOOLS package (Hartman et al. 2008 1 , 


using S2-S5 to reconstruct the light curves of the SN and S 
Due to the fact that S2-S5 are fainter, this method introduced 
noise at the level of 0.01 - 0.02 mags to the light curves of the 
brighter stars, however, we found that the variability signal re¬ 
mains, although it is still affected by systematics. The two lower 
panels of Figures |3}|^present the detrended light curves and the 
differential curve S1-S2, respectively. The latter illustrates sys¬ 
tematics present in the photometry, which are discussed in the 
next Section. Features that are present in all light curves (more 
clearly illustrated in the binned curves) and that correspond to a 
featureless S1-S2 curve provide evidence for rapid variability in 
the SN. 

We next performed a test using artificial stars to assess the 
variability signal seen in the differential light curves. We inserted 
artificial stars to each of the V-band images taken on N3, which 
displays the largest variability, using the appropriate PSF derived 
for each image with the daophot package in IRAF. We inserted 
3 stars at the magnitude of the SN along the disk of the galaxy 
at locations with both smoother and more complex galaxy back¬ 
ground, 3 stars at the magnitude of S2 at the same distance from 
the galaxy disk, and 3 stars at the magnitude and distance of 
S1. We then performed aperture photometry, as described above, 
and constructed differential light curves. These were found to be 
dominated by Poisson noise and not to show any patterns resem¬ 
bling the varibility signal. In particular, the standard deviation 
(cr) of the points for the ’artificial’ SN in all 3 positions tested 
were 4.5-4.8 mmag, while we found cr = 1.6 mmag for SN-Sl, 
cr = 0.0115 mag for SN-S2 and cr = 0.012 mag for S1-S2. The 
artificial star test thus provides additional support for the vari¬ 
ability signal originating in the SN. 


4. Results 

The densely sampled differential light curves shown in Fig¬ 
ures |2}|^ provide a measure of the intranight decline rate of 
SN 2014J. Table[T]presents the decline rates measured each night 
in V and B based on least-square fits to the differential and de¬ 
trended light curves shown in Figures [3}|^ The error bars rep¬ 
resent the root-mean-square (rms) error of the points to the fit. 
Note, the slopes inferred for SN-Sl are identical for the cali¬ 
brated and instrumental light curves. All values are in agreement 


within errors, except for those derived using SN-Sl on Nl, due 
to the larger gaps in the photometry on that night. Adopting the 
values derived from SN-S2, SN 2014J declines with a rate a of 
0.15 + 0.01, 0.09 + 0.01, 0.07 + 0.02, 0.04 + 0.01 mag day^^ in 
the y-band and 0.19+ 0.01, 0.12 + 0.01, 0.13 + 0.02, 0.06 + 0.02 
mag day“^ in the B-banc0 on Nl^, respectively. We therefore 
find the B-band to fade at a faster rate than the V-band, as ex¬ 
pected, and the slope in each filter to vary from night to night. 
The decrease in the slope, observed in both bands, corresponds to 
the onset of the secondary maximum in the near-IR light curve 
( Kasen|[2006 Jack et al. 20151, which is also seen in the light 


curves of SN 2014J presented by Foley et al. (20141, [Amanullah 
P014|l , Marion et al.| ( |2015| l, [Kawabata et aL] ( |2014] l, and 


et al. 


Asha 


i et al.|(|M4r 


More importantly, we report evidence for rapid variability in 
both the V and B-band light curves of SN 2014J. We find that the 
level of variability varies from night to night and is best traced 
by the SN-S1 curve, due to the accuracy in the photometric mea¬ 
surements achieved for these bright stars. N3 and N4 exhibit the 
largest activity with sinusoidal-like variations of amplitude up 
to 0.05 mag, while variability on other nights is typically at the 
0.02-0.03 mag level. A Fourier analysis of the SN-Sl curve did 
not yield a significant periodicity. 


While the precision of our measurements, based on the Icr 
error bars resulting from the aperture photometry, is at the 1.4 
mmag level for SN-S 1 (and at the 0.6 mmag level for the binned 
curve), we must account for correlated sources of error (red 
noise) to quantify the significance of the variability detection. 
Sources of correlated error include the changing airmass, drift¬ 
ing of the stars across the CCD and other instrumental param¬ 
eters. We follow the procedure outlined by |Pont et ah] ( |2006) l, 
which estimates the amount of correlated noise by calculating 
the dispersion cr^ from the binned residuals (after subtracting a 
best fit model, in our case the decline rate) as a function of N 
points and finds the values for the white {(Jw) and red {(Tr) noise 
that best fit the equation: 


N 


+ cr. 


(3) 


The estimated values of cr„, cr^ and cr^ for each night and filter 
are shown in the last 3 columns of Table [T] using the SN-S 1 
data and A = 30 (~ 1 hour), which was selected to represent 
the longest timescale of the detected variability. Note, that the 
values computed with A = 10, 20, 40 were very similar. If rapid 
variability is present, then the ctn values are overestimates, as all 
points of the light curves were used for the calculation. Given 
the derived cr^ values (3-6 mmag), we find the 0.02-0.05 mag 
variations to be statistically significant. At face value, the 0.05 
mag variation in N3 (at HJD=2456707.46-2456707.48) has a 
8.2 (t significance in V and lOcr in B, and the 0.04 mag variation 
in N4 (at HJD=2456708.48-2456708.52) is a l.lcr detection in 
F, and a 6.3cr detection in B. In N2, the 0.02 mag variation at 
HJD=2456706.32 is significant at the 3.8cr level in V and 5.6cr 
in B, while in Nl, the 0.03 mag variation (at HJD=2456705.47- 
2456705.51) is at the 6.2cr (V) and 9.7cr (B) level. The timescale 
of sinusoidal-like variations ranges from 15-60 min. On N3, the 
B-band is found to precede the V-band by ~ 10 min. The B-V 
color gradually increases within each night, displaying an rms 
scatter of about 0.01 mag. 


^ Note, no extinction or atmospheric corrections were applied. 
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Fig. 3; Instrumental B-band (blue triangles) and V-band (green circles) light curves of SN 2014J obtained on N3. The first four 
panels show the differential light curves of the SN with respect to S1-S3 and S5, respectively, followed by a panel showing the 
reconstructed light curves using the trend htting algorithm (TFA) and the differential curve S1-S2. Solid lines correspond to least- 
square hts to the data, while black and red squares correspond to the binned V and B-band light curve, respectively, using a bin 

curves are offset for clarity by the amount indicated in each panel. Dotted lines indicate regions of 

signihcant variabtlity. ^ 
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Fig. 4: Same as Figure]^ but for Nl. On this night all 5 comparison stars were available, therefore the first hve panels show the 
differential light curves with respect to S1-S5, respectively, followed by a panel showing the reconstructed light curves using the 
trend htting algorithm and the differential curve S1-S2. 
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Fig. 5; Same as Figurej^ but for N2. 
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Fig. 6; Same as Figurej^ but for N4. 
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Table 1; Decline rates a (mag day ') and values of o-„, cTr, ctn (mag) based on SN-S1 of the V and B-band light curves of SN 
2014J. 


Night 

Eilter 

o:sN-Sl 

o:sN-S2 

0:sN-S3 

o:sN-S4 

asN-ss 

^SNtfa 

CTw 

(Tr 

cr N 

1 

V 

0.06 ± 0.01 

0.15 + 0.01 

0.17 + 0.02 

0.15 + 0.03 

0.16 + 0.02 

0.14 + 0.02 

0.0086 

0.0045 

0.0048 

2 

V 

0.09 ± 0.01 

0.09 + 0.01 

0.10 + 0.02 

0.08 + 0.03 

0.07 + 0.03 

0.07 + 0.01 

0.0107 

0.0050 

0.0053 

3 

V 

0.10 + 0.01 

0.07 + 0.02 

0.04 + 0.03 


0.08 + 0.03 

0.09 + 0.04 

0.0177 

0.0052 

0.0061 

4 

V 

0.03 + 0.01 

0.04 + 0.01 

0.05 + 0.02 

0.02 + 0.03 

0.03 + 0.02 

0.04 + 0.01 

0.0131 

0.0046 

0.0052 

1 

B 

0.09 ± 0.01 

0.19 + 0.01 

0.24 + 0.03 

0.21+0.03 

0.17 + 0.05 

0.22 + 0.05 

0.0104 

0.0025 

0.0031 

2 

B 

0.14 + 0.01 

0.12 + 0.01 

0.12 + 0.03 

0.10 + 0.03 

0.07 + 0.04 

0.09 + 0.01 

0.0074 

0.0033 

0.0036 

3 

B 

0.16 + 0.01 

0.13 + 0.02 

0.11+0.03 


0.20 + 0.04 

0.15 + 0.03 

0.0186 

0.0036 

0.0050 

4 

B 

0.07 + 0.01 

0.06 + 0.02 

0.11+0.03 

0.09 + 0.04 

0.12 + 0.04 

0.06 + 0.01 

0.0128 

0.0059 

0.0063 


5. Discussion and Conciusions 


The high-cadence, high-precision photometry obtained with the 
2.3m Aristarchos telescope was cmcial in providing evidence for 
rapid variability in the optical light curve of SN 2014J. A ca¬ 
dence of 30 min, such as that of the Kepler Mission, would not 
have been sufficient to detect the sinusoidal-like variations nor 
any variability pattern seen with the 2 min cadence of our obser¬ 
vations. At face value, our photometric precision of 3-6 mmag 
(based on the estimation of red noise) or 2 mmag (based on arti¬ 
ficial star tests on N3) for SN-S 1 allowed for statistically signih- 
cant detections of variability at the level of 0.02-0.05 mag on all 
nights, peaking on N3, implying that the phenomenon is com¬ 
mon or perhaps even ubiquitous in SNe, unless it is produced by 
a mechanism related exclusively to the onset of the secondary 
maximum. 


Theoretical models of SN light curves (e.g. Kasen et al. 
2006[ Sim et al.|2013[ Fink et al.|2014| l do not make predictions 


on such short time scales and therefore cannot be compared with 
the photometry presented here. We propose one or a combina¬ 
tion of the following scenarios for the origin of the variability: 

(i) clumping of the ejecta, possibly caused by structures of in¬ 
termediate mass elements in the outer layers ( Hole et al.|2010| l, 

(ii) interaction of the ejecta with circumstellar material, inferred 
by |Foley et al.| ( |20141 l, (hi) asymmetry of the ejecta, as the ex¬ 
plosion is not expected to be spherically symmetric (as inferred 


from spectropolarimetry, e.g. see review by Wang & Wheeler 
|2008| l, (iv) the onset of the secondary maximum, which corre¬ 
sponds to a sudden decrease in the flux mean opacity due to the 
transition from doubly to singly ionized iron group elements (see 
e.g. |Pinto & Eastman|2000] l. 

Using the photospheric velocity ~2 wee ks after maximum 
(-10,600 km s“*; estimated from Fig. 15 of Foley et al.|2014| l, 
we estimate the SN to have a radius of 104 AU on day 17. The 
0.02-0.05 mag variability corresponds to a fractional change of 
~2-5% in magnitude and flux, respectively, or a ~ 1-2.5% frac¬ 
tional change in radius assuming constant luminosity and that 
the variation is due to a change of radius, i.e. 1-2.6 AU. The ob¬ 
served fluctuation, however, is an average variation of flux over 
the projected surface of the young remnant; the 1-2.6 AU ra¬ 
dius change gives an estimate of the surface area fluctuating. 
Siverd et al.|(2015]) reported relatively high-cadence photometry 


of SN 2014J with the 4.2cm Kilodegree Extremely Little Tele¬ 
scope North (KELT-N), thereby placing a 4.5% (3 (t) upper limit 
on short timescale variations. The evidence for variability at the 
level of 2-5% presented in this work is therefore consistent with 
the results of the measurement with KELT-N. 


In conclusion, we strongly encourage the community to un¬ 
dertake high-cadence and high-precision monitoring campaigns 
of future, nearby and bright supernovae to confirm the presence 
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of rapid variability, determine whether it occurs in both SNe la 
and II light curves, and differentiate between the scenarios for 
its origin. If variability is due to asphericity in the explosion, 
it will provide evidence for distinguishing the nature of SNe la 
progenitors (see Livio & Pringle|201 1 1 . Eurthermore, light curve 
variability, if shown to be ubiquitous, is likely to contribute to 
the scatter in the mean magnitude of SN la calibrators, which 
remains one of the largest factors in the uncertainty of the Hub¬ 
ble constant (e.g. |Riess et al.|201 l| l. In any case, we expect rapid 
variability to provide a new window into the physics of super- 
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Table 2: Calibrated differential photometry of SN 2014J compared to SI in the V and B-bands. 


HJD-2450000 

Filter 

SN-Sl 

CTsN-SI 

6705.30817 

V 

1.231 

0.001 

6705.31330 

V 

1.224 

0.001 

6705.31679 

V 

1.212 

0.001 

6705.31766 

V 

1.220 

0.001 

6705.31853 

V 

1.219 

0.001 



